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ABSTRACT 

We have taken K-band spectra covering 7 cooling flow clusters. The spectra show many 
of the 1-OS transitions of molecular Hydrogen, as well as some of the higher vibrational 
transitions, and some lines of ionized Hydrogen. The line ratios allow us to conclude 
that the rotational states of the first excited vibrational state are in approximate 
LTE, so that densities above 10 5 cm' 3 are likely, but there is evidence that the higher 
vibrational states are not in LTE. The lack of pressure balance between the molecular 
gas and the ionized components emphasizes the need for dynamic models of the gas. 
The ratios of the ionized to molecular lines are relatively constant but lower than from 
starburst regions, indicating that alternative heating mechanisms are necessary. 
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1 INTRODUCTION 

Cooling flows result when the hot, X-ray emitting, gas in a 
cluster is compressed by the gravitational force of the central 
galaxy to a density where it can cool in less than the Hubble 
, time (c.f. Fabian 1994). In many cases 100-1000 M Q /yr con- 
dense out of the hot gas. They present at least two major 
unsolved problems: the ultimate fate of the gas after it has 
cooled, and the source of the copious optical line emission 
in the centre of the flows (Heckman et al. 1989, Voit and 
Donahue, 1997: VD97). 

A new aspect of both these problems was revealed by 
the detection of infrared line emission of Hb molecules from 
the centres of several CF clusters (Mouri 1994, Elston and 
Maloney 1994, Jaffe and Bremer 1997 (Paper I), Falcke et al. 
1998, Genzel et al. 1998). This emission arises from dense 
gas at ~ 2000 K and is not found in similarly radio-loud 
AGNs not associated with CFs. Thus it represents prima 
facie evidence for "cool" material in cooling flows. Its spa- 
tial extent and energetics suggest an intimate association 
with the optical line emitting gas (Paper I). High surface 
brightness molecular line emission is generally concentrated 
within the giant central cluster galaxies (typically in the in- 
ner 10 kpc) in a similar manner to the emission from ionised 
gas. Both can be considered galactic-scale phenomena, even 
though the presence of the gas is statistically correlated with 
the presence of a cooling flow on larger scales (Paper I). 



We do not yet understand the heating mechanism of 
the molecular gas, and cannot accurately estimate the mass 
associated with it. The mass in gas at 2000 K is only about 
10 5-6 Mq but the cooling time from this temperature is 
very short, about one year. Without a consistent thermal 
model, we do not know whether this represents a minor 
mass component (the equivalent of only 1000 years of CF 
mass accumulation) that is kept continually warm, or the 
tail of a cooler but much more massive molecular dog that 
is heated to IR temperatures by shocks, X-rays, or stellar 
photoionization. 

In this paper we describe new K-band IR spectroscopy 
of six cooling flows taken with the United Kingdom InfraRed 
Telescope (UKIRT) on 2 and 3 January, 1997. By studying 
the ratios between the various lines of H2 and between the 
H2 and HII lines we hope to shed light on the thermal state 
of the gas, and ultimately on the heating mechanism. The 
clusters studied were Abell 0335, 478, 1795, and 2029, Hydra 
A, NGC 1275, and PKS 0745-19. They were observed with 
the CGS4 detector in the following configuration: 75 1/mm 
grating; 256(spectral)x 84(spatial) InSb array; slit size:2", 
spatial scale 1" /pixel, dispersion: .0026^/pixel. 

2 OBSERVATION AND REDUCTION 

The observations followed standard IR spectroscopic tech- 
niques. Each obervation consists of a series of 60s frames, 



consisting of four subframes, each shifted 1/2 pixel in the 
spectral direction. The four subframes are combined into a 
single 514x84 frame. 

We calibrated the atmospheric transmission by observ- 
ing, with the same technique, standard stars about every 90 
minutes. We chose the stars from the Carlsberg Meridian 
Catalog (Fabricius, 1993) to be of early F type and within 
10° of the target galaxies. 

For any 1-dimensional target spectrum we used the 
telluric calibrators to subtract the smooth underlying stel- 
lar continuum leaving only non-atmospheric absorption and 
emission features. Our procedure is to find emission lines 
in these spectra, and then correct the fluxes for the atmo- 
spheric absorption, rather than to first correct the spectra 
and then fit lines. For narrow lines, the former procedure 
is less sensitive to spurious lines created by minor atmo- 
spheric fluctuations. These steps are illustrated for the clus- 
ter PKS0745-11 in Figure 1. 

We applied this procedure to the summed 2-dimcnsional 
spectra for each cluster as follows: A continuum model, with 
atmospheric corrections, was subtracted from each row of 
the spectrum. This yields a 2-dimensional spectrum with 
emission/absorption lines only. We then summed the spatial 
profiles of the strongest hydrogen emission lines. The result- 
ing profile was compared to the spatial profile of temporally 
adjacent standard stars to determine the spatial extent of 
the emitting regions. Lastly the rows of the 2-dimensional 
emission line spectrum were weighted by this spatial profile 
and summed to yield a single profile with the optimum sig- 
nal/noise for detecting weak lines, assuming that all lines 
have the same profiles. 

From the uncorrected spectra we determined, by gaus- 
sian fitting, the line widths of the strongest Eb lines. We 
then searched for weaker lines by fitting gaussians of this 
given width at the given wavelengths of a large number of 
H2, HII, and metal lines. Because the position and width of 
the gaussians were fixed, the fitting procedure should have 
no bias toward finding positive peaks. We estimated the un- 
certainties in the fitted fluxes empirically by applying the 
same procedure at 100 random points along the spectrum 
and producing a smoothed r.m.s. average of the fluxes so de- 
termined as a function of wavelength. The method method 
includes residual atmospheric features and usually indicates 
a larger uncertainty in regions near strong telluric lines. 

As a last step, the line strengths are converted to physi- 
cal units by correcting for atmospheric absorption and using 
the standard star magnitudes. In this process an estimate of 
the uncertainty in the absorption is added in quadrature to 
the uncertainty in the gaussian estimate. The resulting flux 
and error estimates are given in Table 1 for each cluster 
except Abell 2029 where no lines were detected during 12 
min of observation. The error estimates do not include the 
uncertainty of the absolute flux scale. This last uncertainty 
is relatively large because few of the Carlsberg calibrators 
have known K-band magnitudes, so their magnitudes were 
extrapolated to K-band from known B and R fluxes using a 
standard F-star spectrum. We estimate the additional un- 
certainty in the absolute line flux scale to be of order 30%, 
but the relative line strengths within one spectrum, or in one 
cluster relative to another, are accurate at the levels quoted 
in Table 1. 



90 
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Figure 2. Calibrated and continuum subtracted spectra for the 
seven clusters studied. Wavelengths are in fi in the rest frame and 
the vertical scales are in 10 — 16 erg cm" 2 s _1 A~ x 



3 REDUCED SPECTRA AND DERIVED 
QUANTITES 

In Figure 2 we present the continuum subtracted and ab- 
sorption corrected spectra for all clusters in physical units 
and in the rest wavelength frame. In Table 2 we present for 
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Table 1. Line fluxes from the studied clusters. Fluxes 


arc in units of 10 16 erg cm 


-2 s- 1 . 




Cluster 


PKS0745 


A0335 


NGC1275 


Abell 1795 


Abell 478 


Hydra A 


Exposure Time (min) 


60 


54 


12 


30 


100 


56 


z 


0.1028 


0.0349 


0.0017 


0.0616 


0.0860 


0.0540 


1-0S(0) (2.2235^) 


3.4 ±5.9 


6.4 ±1.6 


102.0 ±6.0 


-0.6 ±1.4 


2.4 ±1.4 


3.7 ±1.3 


1-0S(1) (2.1218/u) 


24.4 ±0.7 


22.1 ±1.1 


384.3 ± 3.8 


15.8 ±1.0 


13.3 ±0.8 


9.7 ±0.7 


1-0S(2) (2.0338a*) 


9.2 ±0.6 


9.9 ±0.8 


115.1 ±6.2 


7.4 ±0.6 


5.9 ±0.6 


3.1 ±0.6 


1-0S(3) (1.9576/u) 


22.8±0.5 


25.7 ±1.9 


332.3 ± 14. 


17.4 ±0.4 


12.0 ±0.4 


15.5 ±1.6 


1-0S(4) (1.1.8920/*) 


3.2 ±0.5 


6.4 ±11. 


59.2 ± 28. 


3.2 ±8.4 


-1.4 ±0.8 


2.9 ±4.3 


1-0S(5) (1.8358/*) 


9.5 ±0.8 


N.A. 


N.A. 


9.1 ±5.2 


6.9 ±1.3 


6.2 ±7.7 


2-lS(l) (2.2477/*) 


7.2 ±6.0 


-6.3 ±3.0 


34.7 ±6. 


-0.9 ±1.0 


3.3 ±4.3 


3.2 ±2.1 


2-lS(2) (2.1542/*) 


1.4 ±0.8 


-0.1 ±1.0 


10.9 ±4.9 


3.8 ±1.2 


-1.6 ±0.7 


1.3 ±0.9 


2-lS(3) (2.0735/*) 


3.1 ±0.6 


-0.4 ±0.8 


13.1 ±3.6 


2.6 ±0.9 


0.5 ±0.7 


-2.8 ±1.0 


Pa a (1.8756/*) 


36.3 ±0.7 


5.0 ±12. 


300 1 


14.0 ±3.1 


8.4 ±0.6 


17.6 ±1.3 


Br 7 (2.1661^) 


3.7 ±0.8 


-0.6 ±1.0 


47.8 ±4.7 


0.4 ±1.3 


2.0 ±0.9 


1.3 ±1.0 


Br <5 (1.9451/t) 


1.7 ±0.5 


-2. ± 12. 


40.1 ±7.7 


-1.7 ±1.0 


-0.8 ±0.4 


-1.8 ±1.9 



1 The Pao line in NGC 1275 is clearly visible in Fig. 
absorption. The true line flux is highly uncertain. 



2, but lies in a region of strong and variable atmospheric 



each cluster the measured velocity dispersion of the bright- 
est H2 lines, the radial extent of the emission (standard de- 
viation of gaussian fit along slit, deconvolved with stellar 
profile), the total luminosity in the detected H2 lines the 
luminosity in Pace, luminosity in Ha (Heckman et al., 1989; 
White et al. 1994 for Abell 478) and the excitation tem- 
perature T5.3, determined from the ratio of the 1-0S(3) and 
1-0S(1) lines (c.f. Section 4.2). The instrumental dispersion, 
measured on arc and sky lines was 240±10 km s _1 . This may 
be somewhat larger than the effective dispersion for emission 
that does not uniformly fill the slit. In comparing our lumi- 
nosities with those of Heckman et al., determined by narrow 
band imaging, our luminosities should be corrected for light 
blocked by the slit jaws. This can be crudely estimated by 
assuming circular symmetry and calculating the light lost in 
this case from the emission profile along the slit. This leads 
to corrections ranging from 1.6 to 2, but these will be incor- 
rect if the source is highly asymmetric, or if we have missed a 
very large, low brightness component. The latter possibility 
is particularly likely in the nearby cluster NGC 1275. 

The velocity dispersions, corrected for instrumental 
widening, run from 100 to 300 km s _ . These are compa- 
rable to the dispersions measured by Heckman et al. (1989) 
in Ha and are much smaller than the typical velocities in 
the clusters, the stars in the envelopes of cD galaxies or gas 
in the narrow or broad line regions around active nuclei. 



4 DISCUSSION 

4.1 Optical Extinction 

The ratio of Ha to Paa fluxes should measure the extinc- 
tion of the optical light. This ratio, for Case B HII regions, 
should be about 7 (Osterbrock, 1989), and values much be- 
low this would indicate dust around the emission regions. 
Unfortunately the fluxes are measured with different tech- 
niques: slit spectrocopy for Paa and narrow band imaging 
for Ha. Hence the ratio is subject to uncertainties in the 
absolute calibration, and the flux lost outside of the spec- 
troscope slit. For the four high redshift clusters with good 
Paa fluxes, we have taken the ratio LH a /2Lp aa , where the 
factor of 2 corrects for slit losses. The values of the ratio arc 



10, 22, 7, and 4, with an average of 11. While the variation 
confirms the difficulties with calibration, we conclude that 
there is no evidence for strong extinction of the Ha fluxes 
(e.g. more than 1 magnitude). 

4.2 LTE modelling of H2 line ratios 

If the warm gas is sufficiently dense that collisional transi- 
tions dominate over radiative transitions from the relevant 
states, then we can regard the gas to be in local thermal 
equilibrium (LTE). We regard this case as a reference model 
and ask if the measured line ratios are consistent with this 
model. In Table 2, we have given the excitation tempera- 
tures derived from the ratio of 1-0S(1) and 1-0S(3) lines; in 
Figures 3 we plot for each cluster log(F/ag) versus T u . Here 
F is the measured flux in each line, and a, g and T u are the 
Einstein coefficient, statistical weight and energy (expressed 
as a temperture, T u — E u /k) of the upper state of the tran- 
sition. In LTE, of course, log(F/ag) — Constant — T u /T. 

From Table 2 and Figure 3 we conclude that the LTE 
model is a reasonable approximation for the 1-0 lines in all 
clusters, with excitation temperatures of ~ 2000 K. There 
are some lines which deviate significantly more than the for- 
mal errors and these we have investigated individually. In 
some cases we have downgraded the significance of the re- 
sult, because the line appears in a spectral region with a 
poor baseline. This is unfortunately true for many of the 
high excitation, but quite weak, 2-1 lines. By this standard, 
the remaining significant deviations from LTE are as follows: 

For NGC 1275, the 2-lS(l) lies a factor of 3 above the LTE 
model for 1617 K. This is the cluster with the highest sig- 
nal/noise ratios. Krabbe et al. (2000) discuss deeper obser- 
vations of this object, where 2-lS(l),S(3) and 3-2S(3) are 
clearly detected. They find the higher lines to require higher 
excitation temperatures: 2600 K versus 1500K for the 1-0 
lines) and attribute this largely to fluorescent excitation by 
Lya photons. They conclude that for most transitions colli- 
sional processes dominate, requiring n e > 10 5 cm~ 3 . 

For A478 and PKS0745, and marginally for N1275 the 1- 
0S(4) line lies about a factor of 2 below the model. 



Table 2. Quantities derived from analysis of the spectra. 



Cluster 
PKS0745 
A0335 
NGC1275 
Abell 1795 
Abell 478 
Hydra A 



r(kpc) 
1.5 
0.5 1 
0.13 2 
1.1 
1.8 
1.2 



<Tv (km s 1 ) 
320 

370 
300 
280 
260 
370 



Lho 



(erg s x ) 

1.7 10 41 
1.5 10 40 

6.1 10 38 

4.2 10 40 

5.8 10 40 



-Lpaa 
7.3 10 40 

N.A. 

N.A. 
1.0 10 40 
1.2 10 40 



2.3 10 40 9.5 10 39 



1.3 10 42 
5 10 41 

2.1 10 41 

4.4 10 41 
1.7 10 41 

7.2 10 40 



T 5;3 
1771 ± 200 
2446 ± 900 
1628 ± 100 
2217 ± 350 
1693 ± 362 
3 1480±300 



1 A0335-gaussian fit to spatial fit is poor. There is an unresolved core with 80 ± 10% of the flux. The 
remaining 20% extends to at least 3 kpc and has a gaussian fit radius of about 2 kpc. 2 NGC 1275- 
gaussian fit is poor. Unresolved core with ~ 85% of flux. Remainder extends ~ 10" or 3 kpc in 
radius. 3 see text at end of section 4.2 
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Figure 3. LTE Diagrams for the six clusters with detected lines. 
For each cluster the plot shows log F/ag vs. T u with the straight 
line through the 1-0S(1) and 1-0S(3) points. For undetected lines, 
2<r upper limits are given. 



For Hydra A the ratio of 1-0S(3) to 1-0S(1) is so high (1.9) 
that no plausible excitation temperature can be assigned. 
This is probably due to the fact that the 1-0S(3) line falls in a 
telluric absorption band at 2.06/i and the correction for this 
is inaccurate. Excluding this line, the estimated excitation 
temperature for all the other 1-0S lines is 1480 ± 50 K 



4.3 Physical Implications of the H2 emission 

We conclude from the above analysis that for the clus- 
ters studied the rotational states of the first excited vi- 
brational state are close to LTE. This implies a gas den- 
sity high enough that collisional transitions state domi- 
nate over radiative transitions. Using collision cross sections 
from Allison & Dalgarno (1967) and Sternberg & Dalgarno 
(1989), modified using more recent results by Mandy & Mar- 
tin (1993), at T — 2000 K a molecular hydrogen density 
n > 10 6 cnU 3 is required for a 84% thermal population of 
the v = 1, J = 7 level which decays through the H2 V = 1 — 
S(5) line detected in A478 and PKS0745. For the remaining 
clusters, the highest line detected is H2 1 — S(3) which is 
83% thermalized at n — 3 • 10 5 cnU 3 . Using collision cross 
sections from Le Bourlot et al (1999) would even increase 
these densities somewhat. On the other hand, if there is a 
significant electron abundance in the molecular gas, colli- 
sions with electrons would contribute to the excitation, de- 
creasing the required densities. With a fractional ionization 
of 1%, the required densities would be decreased by a factor 
of two, for a fractional ionization of 10% the densities would 
be decreased by a factor of 10. Based on these considera- 
tions it is safe to conclude that the densities of the emitting 
regions are at least 10 5 cm -3 , and possibly even higher. At 
these densities, the line ratios provide little information on 
the gas heating source, since the signature of non-thermal 
excitation processes such as UV-pumping or X-ray excita- 
tion is quenched, except in lines of very high critical density 
which are faint and not detected in the present data. 

The high density and temperature of the gas emitting 
the near-IR lines imply high thermal pressures (nT ~ 10 8-9 
cm _3 K), exceeding by of the order 100 to 1000 that of the X- 
ray gas in the region near the galaxy (T ~ 10 7 , n e ~ 10 _1 ) 
or the ionized optical line emitting gas (T ~ 10 4 , n e ~ 10 2 ) 
as estimated from the [SII] doublet ratio (Heckman et al., 
1989). The latter two phases have comparable pressures. Be- 
cause of this, most models of the optical line emission in the 
literature assume pressure equilibrium with the X-ray phase. 
Some theoretical work has been done on the expected prop- 
erties of the molecular phase in cooling flows (e.g. Fabian 
1994). Our results make it evident that a dynamic, non- 
isobaric model must be used to self-consistently explain the 
ionized and molecular line emission. 

One possible model suggested by the obvious pressure 
imbalance, and the strong linkage of the optical and K-band 
fluxes, is an ablative or "rocket" model of the gas system. 
X- or UV-radiation could heat the surface layers of a cold, 
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gravitaionally bound, molecular cloud to 2000 K. At this 
temperture the gas is no longer gravitationally bound and 
will expand into the surrounding space. As the density de- 
creases the material will ionize and emit the observed optical 
line emission. We will explore this model, and test the pre- 
dicted line ratios, in a future paper. 

4.4 Relation of HII to H2 emission 

In general, the H2 and HII emission in cooling flows are sim- 
ilarly distributed. As noted in Paper I, the H2 and HII emis- 
sion lines correlate well spatially, energetically, and kincmat- 
ically. An extreme example of this is Abell 2029, which to 
date has shown no detectable Ha or H2 IR lines, although it 
is a strong cooling flow and strong radio source. In the two 
nearest clusters, NGC1275 and A0335, the IR line emission 
is dominated by an unresolved, presumably nuclear compo- 
nent. This is partly an artifact of their large angular size 
which excludes some of the extended emission from the slit. 
Both of these cluster have in fact, an extended component, 
as do the spectra of all of the other objects. 

More recently, two of the objects in our sample were 
imaged with H 2 lines at high resolution with HST/NICMOS 
by Donahue et al., (2000). In the images, the H2 emission 
traces Ha emission. This implies that the H2 and HII gas 
are cospatial on scales of a few hundred parsecs or less and 
strongly suggests that the excitation mechanisms of the two 
gas phases are either the same, or closely coupled. 

For the four clusters in Table 2 for which good Paa 
fluxes are available, the ratios of the ratios of Lh 2 / Lpaa 
range from 2.3 to 4.8 (mean=3.4). The single line ratio 
Lpaa/Li-os(i) ranges from 0.6 to 1.8 (mean 1.0; rms 0.4). 
Br7 is generally not detected; the estimated ratio to 1-0S(1) 
is less than 0.15 in all cases. 

The ratios of the HII lines to the H2 lines are lower 
than those seen in "active" comparison galaxies: Seyferts, 
starburst galaxies and Ultraluminous InfraRed galaxies 
(ULIRGs). For all these galaxies, the ratio of Br7/1-0S(1) is 
greater or equal than 1 ( Moorwood &: Oliva, 1988, 1990, 
Puxley et al., 1988, Genzel et al. 1995, Goldader et al., 
1995, Thompson 1995, Vanzi et al., 1998, Genzel et al., 1998, 
Schinnerer et al., 1998, Murphy et al., 1999). We discuss the 
implications of this difference for UV-photoionization and 
shock excitation models below. 

From the consistency of the line ratios, the spatial over- 
lap, and the difference of the ratios from those of other ac- 
tive galaxies, we conclude that in most cooling flows the 
excitation of radiation from the ionized and molecular gas 
are causally related. This may not be true for some galaxies 
with exceedingly active nuclei. Cygnus-A for example shows 
a much higher ratio of Paa/1-0S(1) than we find (Wilman 
et al., 2000). 

4.5 Excitation by stellar UV radiation 

Some authors (e.g. VD97) have suggested that the optical 
spectra of cooling flows can be explained largely by stellar 
photoionization. This cannot be the case for the molecular 
gas detected here. The faintness of Br7 in the cooling flows, 
compared to starburst galaxies, argues against H2 excitation 
by star formation. Photo-ionization by the hottest stars is 
also known to be insufficient to account for the high elec- 
tron temperatures in the ionized filaments of cooling flows, 



which thus already required an additional source of excita- 
tion (VD97). The addition of the IR spectra to the opti- 
cal spectra then potentially requires three excitation mech- 
anisms to be acting simultaneously, and in concert, to yield 
the uniform IR spectra presented here. 

An additional diagnostic is provided by the ortho/para 
ratio, which is expected to be three (the ratio of the sta- 
tistical weights) in situations where the high-temperature 
equilibrium ortho/para ratio is reached (Burton et al. 1992). 
On the other hand, in regions where the H2 lines are UV- 
excited, a significantly lower ortho/para ratio is expected, 
due to spin-exchange reactions with protons and hydrogen 
atoms (Draine & Bertoldi 1996). Indeed, a low ortho-para 
ratio is observed in the extended H2 vibrational emission in 
the starburst galaxy NGC 253 (Harrison et al. 1998) and in 
regions of the well-studied Galactic photodissociation region 
in the Orion bar (Marconi et al., 1998) and the giant extra- 
galactic HII region NGC 5461 (Puxley et al., 2000). The low 
ortho-para ratio in UV-excited regions has been highlighted 
in the models of Sternberg & Neufeld (1999). 

In all of our targets, the v = 1 S(l) (ortho), S(2) (para) 
and S(3) (ortho) lines are detected at high signal-to- noise 
ratio and there is no indication for a deviation from a or- 
tho/para ratio of three in our sample. Together with the 
above discussion, this argues against stellar UV-excitation. 

4.6 Shock Excitation 

A single shock model, where the various emission systems 
are produced in one gas stream that cools after the shock, 
does not predict the observed pressure differences between 
the HII and H2 regions. A shock in cold gas, just strong 
enough to produce the 2000 K H2 lines, would not produce 
the HII lines. A stronger shock could produce the HII lines, 
but, if the H2 lines come from a cooler, post-shock region, 
the pressure there would be lower than directly behind the 
shock. A third alternative, a shock ionizing gas previously 
warmed to 2000 K, would also produce higher pressure in 
the HII gas than in the H2 gas. 

The consistent and low HII/H2 line ratios severly re- 
strict more complicated shock models. As discussed in paper 
I, only a narrow range of shock velocities, below ~ 30 km - , 
could produce the H2 emission without excessive HII lumi- 
nosity, while the characteristic velocities in the centers of cD 
galaxies are well above this. Van der Werf, et al. (1993) dis- 
cuss a similar situation in the merging galaxy pair NGC 6240 
(see further the section on LINERS below). They propose, 
effectively, that dense molecular clouds collide with a low 
density atomic or molecular ISM, forming a fast ionizing 
shock in the low density medium and a slow, C-type shock 
in the molecular medium. They do not explicity calculate the 
ionized to molecular line ratios in this case, and we find the 
existence of such a low density, low temperature, medium 
unlikely at the centre of the cooling flow. 

Alternatively, organized, rotation-like, motions of the 
cloud ensemble may allow the local cloud-cloud collisional 
velocities to be ~ 30 km _1 while the large scale global 
motions are much higher. In this case the origin of the HII 
emission is unclear. It cannot arise directly in the shocks (be- 
cause of the pressure argument) but might arise indirectly, 
for example in a star forming region downstream. 

While the HII and H2 line systems could conceivably 



arise in different regions from different mechanisms (e.g. Hf f 
emission from a AGN narrow fine region and H2 fines from 
shocks in a circumnuciear disk), the iow variation of the 
fine ratios from ciuster to ciuster argues against unreiated 
emission sources. 

We conciude that our resuits can be expiained by sys- 
tems of shocked cfouds, but oniy under rather speciai con- 
ditions, whose iikeiihood requires further investigation. 

4.7 Comparison to LINER spectra 

The reiative strengths of H2(1-0)S(1) and Br7 in our spec- 
tra are remarkabiy simiiar to those found in Lf NER galaxies 
(e.g., Larkin et al., 1998), perhaps not surprisingly as the op- 
tical spectra of cooling flows also have many similarities with 
LINERS (e.g. VD97). Models for the excitation mechanisms 
in LINERS include central AGN photoionization and ioniza- 
tion by an aging starburst (effectively a mix of HII regions 
and supernova remnants, with a relative deficiency in ioniz- 
ing photons accounting for the faintness of Br7), and depend 
on the details of the near-IR spectrum (Alonso-Herrero et 
al., 2000). What produces the H2 vibrational emission in 
LINERs is not clear, except in the class of LINERs where 
shock-excitation dominates the observed spectrum, e.g., in 
NGC6240, discussed above, which has an optical spectrum 
which is typical for shock-excited LINERs (e.g., Fosbury & 
Wall 1979, Fried & Schulz 1983, Heckman et al. 1987). 



5 CONCLUSIONS 

Our UKIRT observations detect multiple lines in six of the 
seven cooling flux clusters observed. In all cases except Abell 
2029, we detect at least four of the 1-0S series, sometimes 
some of the 2-1S series. Paschen a is detected where the 
redshift allows this, but Br7 and 5 are weak. 

The ratio of molecular to ionized hydrogen line emis- 
sion is much higher than those found in starburst regions or 
Seyfert galaxies, but are similar to those in LINERS. The 
ratios of Hq to Paa do not indicate large absorption of the 
former line by dust. 

The 1-0S lines show near-LTE ratios, indicating den- 
sities well above 10 5 cm -3 . This implies pressures in the 
warm gas that are 2 to 3 orders of magnitude higher than 
the surrounding ionized gas, and rule out static equilibria. 

The line ratio and pressure data argue against gas heat- 
ing by stellar UV emission. X-ray heating in non-isobaric 
conditions remains an option, as do shocks under highly con- 
strained conditions. 
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